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The single-mode spin-boson model exhibits behavior not included in the rotating wave approxi- 
mation (RWA) in the ultra and deep-strong coupling regimes, where counter-rotating contributions 
become important. We introduce a symmetric rotating wave approximation that treats rotating 
and counter-rotating terms equally, preserves the invariances of the Hamiltonian with respect to 
its parameters, and reproduces several qualitative features of the spin-boson spectrum not present 
in the original rotating wave approximation both off-resonance and at deep strong coupling. The 
symmetric rotating wave approximation allows for the treatment of certain ultra and deep-strong 
coupling regimes with similar accuracy and mathematical simplicity as does the RWA in the weak 
coupling regime. Additionally, we symmetrize the generalized form of the rotating wave approxi- 
mation to obtain the same qualitative correspondence with the addition of improved quantitative 
agreement with the exact numerical results. The method is readily extended to higher accuracy if 
needed. Finally, we introduce the two-photon parity operator for the two-photon Rabi Hamiltonian 
and obtain its generalized symmetric rotating wave approximation. The existence of this operator 
reveals a parity symmetry similar to that in the Rabi Hamiltonian as well as another symmetry 
that is unique to the two-photon case, providing insight into the mathematical structure of the two- 
photon spectrum, significantly simplifying the numerics, and revealing some interesting dynamical 
properties. 
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The Hamiltonian of a two-level system coupled lin- 
early to a quantum harmonic oscillator [l|, Q is one of 
the most studied models in quantum mechanics and is 
still unsolved despite over fifty years of nearly continu- 
ous effort. The model has found applications in many 
fields ranging from molecular chemistry to circuit quan- 
tum electrodynamics, which is evident when one lists the 
different pseudonyms for it in the literature: the Jaynes- 
Cummings model (without the rotating wave approxi- 
mation) |3, H| and the single-mode spin-boson (H-0] in 
condensed matter physics, the Rabi Hamiltonian in quan- 
tum optics the molecular dimer-oscillator in chemical 
physics n|, and the spin-oscillator 11 1 and qubit- 
oscillator in quantum information and circuit QED (l2| - 
[T3 |. Many analytical approaches have been developed 
for small and (more recently) large coupling between the 
system and the oscillator, with arguably the most well- 
known being the rotating wave approximation (RWA). 
In this article we introduce and discuss approximations 
using similar techniques as those used in obtaining the 
RWA, but in parameter regions complementing it. 

The RWA Q is designed to work well in the case of 
weak coupling between the two-level system and the os- 
cillator. Additionally, previous approximations [HI, [l6| 
have been linked to a generalized form of the RWA in 
Ref. [l?]]- However, recent experiments [l8[ have mo- 
tivated the study of this model in parameter regions 



which have not been thoroughly explored in the past 
and have shown that the RWA breaks down in those re- 
gions [l3|, [l9| . Specifically, a number of recent theoretical 
studies have shown that contributions of counter-rotating 
terms, which are ignored in the RWA, prove important in 
these parameter regions and exhibit d yna mical behavior 
different from the weak-coupling case 

Addi- 
tionally, counter-rotating terms dominate the short-time 
dynamical behavior for some parameter regions, leading 
to important Zeno and anti-Zeno effects that are not re- 
produced by the RWA [H-Hl]. 



The Hamiltonian of a two-level system coupled 
quadratically to a quantum harmonic oscillator, the two- 
photon Rabi Hamiltonian, has also been studied within 
the RWA H3- Limitations of the RWA have been out- 
lined for this system [28[, but so far limited effort has 
been directed to studying it outside of the RWA [29i — 
|32| . This Hamiltonian arose in quantum optics as a phe- 
nomenological model for a three-level system interacting 
with two photons [28|, [3(3, [33| and is also relevant in mod- 
eling pure dephasing in crystals (34|. As opposed to a 
displacement in position in the case of the Rabi Hamilto- 
nian, the coupling in the two-photon Rabi Hamiltonian is 
through frequency displacement or "squeezing" [35l [36j . 
With bi-exciton effects and two-photon processes occur- 
ring in experimental systems [37j, more work is needed 
to determine whether this Hamiltonian can successfully 
model these effects. 
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The two Hamiltonians (m = 1,2) can be written in the 



form 



where 



H m =Ljtfb + Ja x + Xa z (b rn + (fet) 



(1) 



where H\ is unitarily equivalent to the single-mode spin 
boson/Rabi Hamiltonian and H2 is the two-photon Rabi 
Hamiltonian. Here, <t,; are the usual Pauli matrices [38| . 

and b are the boson raising and lowering operators, u is 
the harmonic oscillator frequency, J is the coupling of the 
two-level system, and A is the coupling strength between 
the two-level system and the harmonic oscillator. The 
recent regions of interest for H± include the ultra-strong 
(A > 0. la;) and deep-strong (A > w) coupling regimes 



In this work, we introduce a symmetric form of the ro- 
tating wave approximation (denoted as S-RWA) that in- 
cludes an equal amount of rotating and counter-rotating 
terms. The S-RWA provides analogous physical insight in 
the off-resonance ultra and deep-strong coupling regions 
to that provided by the RWA at resonance in the weak 
coupling limit. Similar to previous use of the RWA in the 
weak-coupling limit [3j|, one can employ the S-RWA to 
extract important features of the dynamics in the strong 
coupling limit. In Section HI we discuss a mathematical 
feature of Hi conserved by the S-RWA: the invariance of 
its energies under change of sign of the coupling param- 
eters J and A. In Section [Til wc symmetrize the RWA 
and describe the parameter regions in which it is most 
applicable. In Section [ill! we provide a generalization of 
the S-RWA. Similar to the generalization of the RWA to 
larger coupling [13], the analogous generalization of the 
S-RWA (denoted as S-GRWA) extends it to regions of 
small coupling (sketched in Table [TJ) . The S-GRWA can 
also be extended to higher levels of accuracy if needed, 
a feature that is not as directly evident in the general- 
ized RWA. Finally, by introducing the two-photon parity 
operator in Section IIV1 we are able to apply analogous 
methods and obtain an S-GRWA for the two-photon Rabi 
Hamiltonian H-2- We show that, apart from respecting 
the same symmetries present in H\ , the two-photon case 
contains two independent manifolds and the two-photon 
parity operator maintains parity symmetries on each of 
them. 



U r (</>) = exp (i(f>a y /2) . 



(2) 



is equivalent to letting J —> — J. As expected, the set 
of eigenvalues of this Hamiltonian, {e ± J}, is invariant 
under this transformation. The analogous boson problem 
is the displaced harmonic oscillator 



H+ = ujtfb + X(b + 6+) 



(3) 



which can be diagonalized with the unitary displacement 
operator (3f| 



V(X/uj) = exp [(&- 6+) X/uj] , 



(4) 



yielding the set of energies {ujN — ^-}^ =0 - Applying the 
well-known bosonic parity /reflection operator, 

K = exp (inb%) , (5) 

is equivalent to letting A — > — A. In this case, the in- 
variance is more pronounced: while the two-level system 
involves a re-ordering of the index {e ± J} — > {e =p J}, 
the energies of both the original and transformed boson 
Hamiltonians are equal for each N due to their depen- 
dence on A 2 . 

For H m from Eq. ([1]), both of these invariances are 
present and can be confirmed numerically. For the Rabi 
Hamiltonian Hi, applying the transformations 1Z and 
U r (w)Ti. is equivalent to letting A —> — X and J —> — J, 
respectively. It is important to note that the invariance 
with respect to J is not a mere convention: for any real 
J, symmetrically correct sets of energy approximations 
have to be the same for both J and — J. The exact ener- 
gies of Hi could therefore be dependent on even powers of 
J and/or contain a "±" splitting for odd powers. While 
these concepts may seem trivial, the invariance with re- 
spect to J, while maintained for some newer approxima- 
tions HHim, is broken for the RWAJI, GRWA [13], 
and other approximations [fl [ll], EH, Elj . as discussed 
below. 



II. SYMMETRIZING THE ROTATING WAVE 
APPROXIMATION 



I. SYMMETRIES IN BOSONIC AND 
FERMIONIC SYSTEMS 

Changing the basis for a Hamiltonian H is equivalent 
to applying a unitary transformation U to obtain a trans- 
formed Hamiltonian H = U^HU, with the eigenvalues of 
H preserved [4(|. Therefore, if there exists a unitary 
transformation such that applying that transformation 
corresponds to changing the sign of a parameter of H, 
then the set of eigenvalues must be invariant under this 
change of sign. 

We note two relevant examples of this property The 
first is the degenerate two- level system: H = e+Jcr x . Ap- 
plying the two-level fermionic reflection operator U r (%) 



The RWA consists of rotating the system by U r {ir/2) 
and then transforming into the interaction picture with 
respect to the interaction-free Hamiltonian H^ ' = ujb jf b+ 
Jcr x . Assuming J > 0, the terms which rotate at frequen- 
cies of u + 2 J are assumed to oscillate much faster than 
the terms rotating at u — 2 J and are thus ignored. Note 
that the counter-rotating wave approximation (CRWA) 
ignores the slower-oscillating terms, i.e., ui — 2 J, and ob- 
tains the excited state energies of the RWA with the sign 
of J changed [H]. The RWA energies arc E$ WA = - J 
and (for N > 0) 

eT± a = ( n + \) w± \Z(r- J ) 2 + A2 {N + 1) - (6) 
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Note that the RWA energies are proportional to a dif- 
ference between ^to and J, known in quantum optics as 
the detuning frequency. Thus, any given set of RWA (or 
CRWA) energies is not invariant under changing the sign 
of J due to the symmetry breaking that these approxima- 
tions cause. An additional effect of this is the existence 
of an isolated ground state that does not depend on the 
coupling A, which has been pointed out previously [H,|43jj. 
Using the methods outlined below, we obtain the sym- 
metric RWA energies (with p = ± and for N > 0) 



E 



S-RWA 
N,p,± 



2N+ \ ' 



(7) 

The extra "parity index" p is precisely the reason for in- 
variancc under the sign of J. The isolated ground state 
disappears and the S-RWA ground state, Eq~^_ a , is of 
the same accuracy as the excited states. It should also be 
noted that the excited RWA (CRWA) states for even N 
are exactly the positive (negative) parity S-RWA states, 
demonstrating that the symmetrization maintains the 
relative mathematical simplicity of the RWA while also 
including rotating and counter-rotating terms equally. 



A. Obtaining the S-RWA 

To obtain Eq. j7|) , note that Eq. ([T]) can be written ex- 
plicitly in the spin-i basis |±) (with |+)(+| being the top 
left entry) , where each of the four entries are operators on 
the boson Fock space. In this form, both Hamiltonians 
are 



Hn 



H+ J 
J H„, 



(8) 



where H^ = ub^bzLX (b m + (6^)™) arc symmetric matri- 
ces in the Fock (number-state) basis. Instead of rotating 
them by U r (n/2) as prescribed by the RWA, we introduce 
a non-local unitary transformation 



1 



_ ( 1 - V " 

V2 \ K i 



(9) 



that will diagonalize these Hamiltonians in the spin sub- 
space, with the unitary V m determined below. In order 
for the two off-diagonal terms of H m = U} n H m U m to 
vanish, the following two equations must be satisfied: 



h: 



V = 1. 



(10) 

(11) 



The second equation shows that V m is in fact a parity 
operator: V m = V.^ 1 = V) n . If such a V m exists, then 
the transformed Hamiltonian H m is diagonal in the spin 
subspace: 



H m = 











H: 



(12) 



Various forms of this transformation for Hi have been 
proposed independently 0,[i|[2l|,[43|,[44| and the transfor- 
mation is extendable to the multi-mode case. For m = 1, 
Ref. [l[ gives V\ = TZ = exp (iirtfb) and the Hamiltonian 

Hi consists of the diagonal entries 

Hf =ujtfb± A (6 + 6 f ) ± JTZ (13) 

which are both Hamiltonian operators on the boson Fock 
space. The H 1 provide analytical and numerical advan- 
tages over the original form of the problem [|| |^, [U liH - 
WJ\ . Additionally, they differ only by the signs of the 
parameters and correspond to parity-related subspaces 
@ or "parity chains" [2l| . Writing in the Fock state 
basis | A) (where N = 0, 1, 2, ...), we make the approxima- 
tion by truncating them into two sets of 2-by-2 diagonal 
blocks 

2Noj ± J ±W2N + 1 
iAvWTT (2N+1)ojtJ 

The eigenvalues of these blocks are the symmetric RWA 
energies given in Eq. (O, a result that we call the S- 
RWA. The S-RWA is plotted with the RWA in Fig. [Hat 
J = \lo. At a given A in the ultra-strong and (even more 
so) deep strong coupling regimes and in the off-resonance 
cases (J < iw), energies of opposing parity approach 
one another with increasing N and the exact numeri- 
cal spectrum of the spin-boson can be divided into two 
nearly degenerate columns of differing parity. In these 
regimes, and in this limit, parity no longer determines 
the energy but merely defines the column to which a 
given energy belongs. It is thus applicable and appropri- 
ate to analyze each parity case separately at deep strong 
coupling [2l|. The failure of the RWA at ultra-strong 
coupling can be seen in Fig. [1] In the limit of large 
A, the respective upward and downward-sloping RWA 
energies are equally spaced while the S-RWA energies 
of opposing parity pair up and correctly approach one 
another (E%™^ « E%™%). Additionally, the RWA 
ground state is surpassed by the first excited state at 
A ~ 0.7 (black circle), resulting in an un-physical flip- 
ping of the well-defined ground state parity not present 
in the exact results. On the other hand, the S-RWA 
ground state is maintained throughout the entire param- 
eter range shown. The trade-off to the improved behavior 
at large coupling is the presence of un-physical crossings 
of the negative parity S-RWA excited state energies at 
A < 0.5. In order to better understand the origin of these 
features and the physical relation of the S-RWA to the 
RWA, we examine the two approximations analytically 
in the next section. 



B. Comparison with the RWA 

When performing the RWA, Hi is written in the eigen- 
basis of a x coupled to a Fock state: 

\±x,N) = -±=(\+)±\-))®\N). (14) 
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There are only two types of off-diagonal terms in this 
basis: resonant (+x, N\H±\ — x, N + 1) and off-resonant 
(—x,N\Hi\ + x, N + 1) (and their respective conjugate 
transposes). The RWA (CRWA) ignores the off-resonant 
(resonant) contributions and keeps the resonant (non- 
resonant) terms [T3|. In the case of the S-RWA, the 
transformation U\ allows a different partitioning of these 
off-diagonal terms, making it possible to symmetrically 
ignore half of the resonant and half of the off-resonant 
contributions. Rotating by U% is equivalent to writing 
Hi in the basis 



\4>±,n) 



-L(\±,N)±(-lf 



T,N) 



(15) 



The original basis | ± x, N) is thus mixed in a such a way 
that (4>±,n\Hi\^ t>m ) = for all N,M. We can there- 
fore write Hi separately in terms of \<P+.n) and \4>-.n), 
resulting in H^ and H^ , respectively. Similar to the 
RWA, we now observe that the off-diagonal terms can 
be partitioned into two types: (4 > ±,2n\Hi\(/)± i 2N+i) and 
(4>±,2N+i\Hi\(f)± i 2N+2) ■ The first type consists of even- 
to-odd Fock-space transitions which couple |±x, 2N) and 
| =F x, 2N + 1) while the second are odd-to-even transi- 
tions coupling | ± x, 2N + 1) and | =p x, 2N + 2) (up to 
a ± overall phase) . The S-RWA keeps the first type and 
ignores the second, which is equivalent to removing half 
of the resonant and half of the off-resonant terms. 



The RWA is an approximation that is valid for the 
resonance case (w « 2J) and at small coupling since ro- 
tating terms dominate the long-time dynamical behavior 
in this region: ui — 2 J <C ui + 2 J. One the other hand, 
the S-RWA is the analogous approximation in the com- 
plementary off- resonance (uj ^> 2 J) and deep-strong cou- 
pling region, where u — 2J « uj + 2J. In this regime, 
both the rota ting and counter-rotating terms are equally 
as important |24j . This leads to novel short-time dynam- 
ics and the Zeno effect in the single [25[ and multi-mode 
spin-boson systems 0, H^], topics of significant interest. 
While neither the RWA nor the S-RWA energies from 
Eq. ([7]) are particularly accurate in quantitatively re- 
producing the full numerical spectrum, generalizing the 
S-RWA, in a method analogous to generalizing the RWA, 
produces improved agreement with the exact results at 
resonance and at large coupling. 



III. SYMMETRIZING THE GENERALIZED 
RWA 

The generalization of the RWA involves a simple 
change of basis involving T> from Eq. (jj) prior to per- 
forming the analogous 2-by-2 matrix truncation per- 
formed in the RWA. The GRWA thus includes behaviors 
of the adiabatic approximation [l7j , extending the valid- 
ity of the RWA to arbitrarily large couplings. Since the 
S-RWA is valid for only large values of the coupling, the 
symmetric generalized RWA extends the validity of the 
S-RWA to arbitrarily small coupling in the off-resonance 
regime. The regions of applicability of the respective 
extensions of the RWA and S-RWA are qualitatively de- 
picted in Table HI 

The derivation of the S-GRWA is similar to that of the 
GRWA since a similar change of basis is performed in 
both. The diagonal Hi is now transformed with Ux> = 
T> ( — <j z ) , removing the linear bosonic coupling terms in 
Hf, obtaining [45[ 



jtfb - \ 2 /uj ± JTZV (±2A/w) 



(16) 



These Hamiltonians are written solely in terms of the 
number operator and two forms of the displaced parity 
operator [48j ]. which differ only in the displacement di- 
rection. To write them in Fock space, we need the Fock 
space matrix elements of the displacement operator [35| , 
D m ,n = (M\V(2\/u)\N), which artQ 



D 



M,N = 



N-M 

L 



N-M 
M 



4A 2 ' 

T3- ) ■ ( 17 ) 



To obtain the symmetric GRWA, we write Hf in Fock 
space and, mimicking the method utilized in the gener- 
alized RWA [17[ , truncate them to 2-by-2 block-diagonal 

A 2 

form. Excluding the "— — " term, the iVth block is then 
given by 

2NLO ± J D 2 N,2N J D 2 N,2N+l 

J D 2 N,2N+1 (2A + 1) UJ =F J D 2 N+1,2N+1 

"(18) 

Since the transformed matrix is full (as opposed to hav- 
ing many zero entries as is the case of the GRWA) , there 
is the additional liberty of choosing the size of the matrix 
truncation. Surprisingly, energies identical to the 1-by-l 
block diagonal truncation have been obtained previously 
via other methods HE S3- We also note that the 2-by- 
2 and 4-by-4 truncations were performed in a very similar 
treatment Q, but the negative parity excited state en- 
ergies were discarded. The 2-by-2 truncation performed 
above is of the same order of accuracy as the GRWA, 
and diagonalizing Eq. (|18[) gives the following S-GRWA 
energies: 



1 Here, L^f M (x) is an associated Laguerre polynomial, which 
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1 /2 

± ( \ \--v{^ [Lin (4A> 2 ) + L 2N+l (4A 2 /^)] } 2 + u2 ff + 1} e^"' [2^ (4A> 2 )] 2 ) 



As with the positive parity (p = +) S-RWA energies be- 
ing equal to the even RWA excited state energies, the 
positive parity S-GRWA energies are exactly the even 
GRWA excited state energies. The numerical (solid 
black) and approximated energies (dashed) are plotted 
against the coupling A in Fig. [^a) for the resonance case 
2 J = uj = 1, where differences between them are most 
evident. At couplings A > 0.5, the agreement is almost 
identical with the GRWA, with the exception of the im- 
proved ground state of the S-GRWA. The applicability 
of the S-RWA is vastly extended to the regions of smaller 
coupling A < 0.5, but the GRWA still maintains better 
agreement with the exact results in that region. The neg- 
ative parity energies, alternating by two with the positive 
parity energies, are not in good agreement with the ex- 
act results at small coupling. In order to correct this, 
one can easily program the 4-by-4 matrix truncation into 
Mathematica, which has been done in Fig. [2jb) for a 
larger coupling range. Here, half of the artifacts of the 
negative parity disappear, occurring in two out of every 
eight excited state energies as opposed to every four. In 
both truncations, these artifacts decrease with increasing 
energy level and disappear in the off-resonance regime, 
making the 4-by-4 truncation virtually indistinguishable 
from the exact energies in that regime (not shown) . The 
initial physical reasoning for the respective RWA and S- 
RWA approximations suggest that the GRWA should be 
used for the small coupling regime while the S-GRWA, 
with its respect for parity and a corrected ground state, 
should be used for A > 0.5. 



equal to zero: 













±A\/2 












1 1 

1 UJ 1 

i i 









| ±AV6 









2uj 














: ±aV6 ; 









1 1 

i i 



v ••• •••/ 

The square (dashed square) boxes represent the even 
(odd) manifold couplings; all other couplings are zero. 
This is equivalent to stating that (2iV|F^|2M + 1) and 
(27V + \\Hf\2M) are zero for all N,M. Thus, one can 
reorder the basis and write H 2 first in the even basis 
and then in the odd, converting the pentagonal into two 
tridiagonal matrices: 

/ ±\V2 \ 
±A\/2 2oj '■■ 

o '•. '■■ 

uj ±\V& o 

±A\/6 3w '•• 

V o -J 



IV. TWO-PHOTON RABI HAMILTONIAN 

A. Even-odd separation 

Consider the two-photon case described by Hi from 
Eq. (fTJ). We first note a symmetry of this system that 
is not present in Hi, namely that Hi does not couple 
even and odd Fock states. This has been discussed in a 
different context [32| and stems from the nature of the 
quadratic coupling in the squeezed harmonic oscillator 
H 2 from Eq. ([5]). Writing H^ in the Fock state basis, we 
arrive at a pentagonal matrix with two middle diagonals 



can be denned for all values of integers M and N |35| . Since 
V(-2\/u) = I>t(2A/uj), the elements (M\V(-2\/uj)\N) = 



It is important to stress that the even and odd mani- 
folds are orthogonal and independent of each-other. This 
property admits an interesting observation: if the initial 
wavefunction for this system contains only even (odd) 
Fock states, then it's evolved wavefunction will remain 
in the even (odd) Fock state manifold. Additionally, if 
the initial wavefunction is comprised of both even and 
odd Fock states, then those states will evolve separately 
and will not exchange populations with one another at 
any time. 

This property of Hi is conserved under unitary trans- 
formation, is independent of the values of the parameters, 
and will be used to obtain symmetric GRWA energies for 
this system. As a result, there will be one set of ener- 
gies for each manifold, as opposed to just one set for Hi . 
In addition to this even-odd decoupling, the invariances 
with respect to J and A that hold for Hi also hold for 
Hi, as shown below. 
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B. Two-photon parity 



We introduce the two-photon parity operator: 



T = exp 



iir-tfb(tfb-l) 



(20) 



Acting on a Fock state \N), the operator can be written 
as follows: 



N-l 



T\N) = (-1) 



,N(N-1) 



\N) =exp hV^i \N). (21) 



j=o 



Note that the eigenvalue of T for a state |JV) is a product 
of the parities of the previous N— 1 states. Since ^N(N— 
1) is always an integer, it is easy to prove that T = = 
T~ x ■ It is evident that T maintains the same behavior 
for the manifold of even and odd Fock states as 1Z does 
for the set of all states: 

oo 

T = ('if (\2N)(2N\ + \2N + 1)(2N + 1|) . (22) 

N=0 

This is not accidental, given that H 2 is decoupled within 
the even and odd state subspaces as mentioned above. 
Thus, there is a parity on each of the even/odd subspaces, 
giving a total of four parities: {— l e , — 1 0) +le, +l o }0 

The two-photon parity T commutes with 1Z and its 
actions on the raising and lowering operators are: 



Tb = -bllT 



(23) 



We know from Eq. CGI that Wall = -b and TZb^TZ = 
—tf. Thus, T satisfies Eq. (fTU)) and anti-commutes with 
the quadratic raising and lowering operators: 



Tb 2 = -b 2 T 
Tb r2 = -tf 2 T. 



(24) 



The two-photon parity T thus satisfies both conditions 
imposed by Eqs. (fT0|) and (|TT|) . implying that it is indeed 
the parity operator V2 needed to diagonalize the Hamil- 
tonian in the spin subspace using the transformation U2 
from Eq. ^ . Similar to Hi , the two-photon Rabi Hamil- 
tonian energies are also invariant under change in sign of 
A and J. Applying the transformations T and U r (ir)T 
is equivalent to letting A — > — A and J — > — J, respec- 
tively In analogy with the displaced harmonic oscillator 
Hi from Eq. ((3j), the squeezed harmonic oscillator H^ 
is diagonalized by the unitary squeeze operator (ioj 



S{a) = exp [\a(b 2 -tf 2 )} 



(25) 



2 We note that the Fourier operator \TR. has been labeled as the 
"pseudo-parity" of H2 l29t |30| . |48H because the true parity has 
previously been unknown. 



producing the set of energies {ujN — ^ (oj ■ 
where 



& = uo\Jl - (2A/w) 2 



The parameter a in S(a) satisfies the conditions 

sinh 2a = 2\/Cj 
cosh 2a = oj/uj. 



(26) 



(27) 



Equation (f2"o]) implies an upper bound for the squeezed 
coupling: A < ^w, and thus < Co < lj. More impor- 
tantly, the energies of the squeezed harmonic oscillator 
are invariant under change of sign of A, as is the case 
of the displaced oscillator. Thus, H2 preserves all of the 
invariances of Hi and contains the additional ability to 
be divided into even and odd Fock state manifolds. 



C. S-GRWA for the two-photon case 

We perform the same steps as with Hi in Section Hill 
except that now the two-photon parity is used (1Z — > T) 
and squeezed as opposed to displaced Fock states are 
implemented (V — > S). First, one obtains Hi = U2H2U2 
with diagonal entries 



Hf = ujtfb ± A (b 2 + 6 t2 ) ± JT. 



(28) 



The Hamiltonian H 2 is thus divided into two parity- 
related subspaces. Furthermore, each of the parity spaces 
can be divided into even and odd subspaces, giving one 
subspace per each of the four parities of T. The four 
operators are tridiagonal matrices in Fock space and sig- 
nificantly decrease the effort required to obtain accurate 
numerical energies for H2, paving the way for the efficient 
application of diagonalization schemes [46j . Additionally, 
this four-fold separation reveals that diagonalizing H2 
is equivalent to determining the spectrum of an infinite 
tridiagonal matrix with all three diagonals approaching 
infinity, as is the case for Hi Q. The differences are 
that the diagonals increase faster and there are four as 
opposed to two matrices due to the extra even/odd man- 
ifold symmetry. Finally, the Hamiltonian Hi commutes 
with the parity operator lZa x , a combination of the inver- 
sion operators for the respective spin a displaced boson 
systems. The Hamiltonian H2 commutes with both Ta x 
and 1Z. The reflection operator determines the manifold 
(e or o) while To x determines the parity on each manifold 
(±1), giving the four parities mentioned above. 

Utilizing 2-by-2 block diagonals of the four tridiago- 
nal matrices in Eq. (f2"5]l gives the symmetric KWA for 
H2 ■ In order to obtain the more accurate generalized ver- 
sion of these energies, we apply Us = S (aa z ) to further 



transform H^ into: 



Hf = Cvtfb - i (u - w) ± JTS (±2a) . 



(29) 



G 



ifolds, and truncating to 2-by-2 diagonal form (although 
higher truncations are again possible) gives the following 



Here, we use the fact that TST = , just like 1ZV7Z = 
TV in the linear case. The squeeze operator is also sep- 
arable into even and odd subspaces, and its Fock space four matrices (excluding the "—5 (w — £>)" term) 
matrix elements Sm,n = (M\S(2a)\N) for real a > 
are0 



S- 



2M,2N = 



(2M)! W M+N \u 



(30) 



52M+1 - 2W+1 ~ V( ( 2Aftl)! (S) P m+^+i(S)- (31) 
Writing in the Fock basis, dividing into the two man- 



ANu! ± J 04JV.4JV JS4NAN+2 

JS4NAN+2 [4N + 2) W =F JS l 4iV+2,4JV+2 

(47V + 1)0 ± JSW+ 1.4JV+1 J54AT+l,4iV+3 

JS4JV+1.4AT+3 (47V + 3) ill =p J54iV+3,4iV+3 



Diagonalizing these matrices gives two sets of symmet- 
ric GRWA energies for the two-photon Rabi Hamiltonian 
(with parity index p = ±): 



J 



E 



iTP ,even 
N,p,± 



^47V + (D - ^ +p^v / ^7^[- p 4JV (w/w) - Puv+2 (w/w)] 



(32) 



wJ 2 



1/2 



TP. odd 



± N W - p-VS/w [iW («/w) + P 4W+2 (tf/w)] J + (4jV + 2) (4JVTT)^ J 

^47V + a) - I + [^4W+i (w/w) - Pyv+3 (w/w)] (33) 

± ( I a) - x/w/w [-P4W+1 (w/cj) + P4N+3 (w/w)] I + 



a) J 2 



(47V + 3) (47V + 2) cj 



These energies are plotted in Fig. [3] against the numer- 
ically obtained curves (black) for (a) J = cj = 1 and 
(b) 2 J = w = 1. Since there are two orthogonal man- 
ifolds (even and odd) and a parity (±) for each, there 
are four different types of curves. Red (blue) denotes the 
even (odd) manifold and large (small) dashes are used 
for positive (negative) parity on each respective manifold. 
Although the numerical energies have the same parities 
and come from the same respective even/odd manifolds 
as the analytical approximation, all numerical results are 
in black for an easier visual comparison to the approxi- 
mation. In the resonance case [Fig. [3] (a)], the respec- 
tive even and odd manifold ground states are modeled 
well, with the ground state of the whole system being 
the ground state of the even manifold, E^ l P j e ^ en , for all 
values of the coupling. However, if one is to have an 
initial state consisting of strictly odd Fock states, then 
that state would evolve in a system whose ground state 



is that of the odd manifold, i.e., E, 



TP, odd 
0, 



. Similar to the 



spin-boson case, the negative parity energies are not very 
accurate at small values of A. However, the errors in en- 
ergies for small coupling decrease with increasing energy 



3 Here, Pjg-i-if(. x ) is an associated Legendre polynomial, which 
can be defined for all values of integers M and N. Since 
S(-a) = S^(a), the elements (M\S(-a)\N) = S N , M - These 
were obtained using Ref. |49|| and Eq. H27I I; they are in agree- 
ment with previous efforts l5Cj - 



I 

level, just like for Hi . Another similarity to Hi is that in- 
tersections occur between states for each respective man- 
ifold and between manifolds. The exact determination of 
the parity of each curve allows one to clarify the relative 
meaningfulncss of intersections between them, building 
on previous efforts which have analytically determined 
these " Juddian points" J2j| [32[ . Additionally, the form 
of the Hamiltonian in Eq. (|2l?|) allows one to see that H2 



approaches 



■iw as A 



In the off-resonance case 



of 2 J = lu = 1 [Fig. [3] (b)], the results are significantly 
better and the different crossing pattern between parities 
and manifolds is reproduced by the analytical energies. 



V. CONCLUSION 

We have used a previously known unitary transforma- 
tion to diagonalize the single-mode spin-boson Hamilto- 
nian in the spin subspace and obtain symmetric versions 
of the regular and generalized rotating wave approxima- 
tions that respect the symmetries of the Hamiltonian and 
include an equal amount of rotating and counter-rotating 
contributions. Additionally, we have devised an anal- 
ogous unitary transformation for the two-photon Rabi 
Hamiltonian and obtained its respective symmetric gen- 
eralized rotating wave energies. The symmetric rotating 
wave approximation (S-RWA) allows the short-time dy- 
namics of the single-mode spin boson and two-photon 



7 



Rabi Hamiltonians to be analyzed in the ultra and deep- 
strong coupling regimes with the relative numerical sim- 
plicity of RWA-type approximations. The truncated S- 
GRWA energies shown in this work and the two other 
proposed higher-order truncations can be used to pro- 
vide a highly accurate quantitative picture of the details 
of the dynamical behavior for both Hamiltonians. Both 
the S-RWA and S-GRWA can be applied to model exper- 
iments where there is a large coupling between the spin 
and bosonic systems and where counter-rotating terms 
and parity are important [H, HIM [H, [H[. 

The invariances with respect to the signs of J and 
A present in the uncoupled Hamiltonians and 
(m = 1,2), respectively, demonstrate that the concept of 
"parity symmetry" for the coupled systems Hi and H2 
stems from the symmetries of these uncoupled Hamilto- 
nians. Bosonic Hamiltonians with m > 2 represent 
anharmonic oscillators and cannot be completely diag- 
onalized by operators similar to V and S [40| • As a 



result, the generalized treatment of the approximation 
in this work cannot be applied to the ill-defined [5l[ m- 
photon systems H m with m > 2. Nevertheless, some 
of the anharmonic boson Hamiltonians do contain 
a definite parity symmetry: 1Z is a parity operator for 
the odd anharmonic cases m = 3, 5, 7, ... and T is a par- 
ity for m = 6, 10, 14, ... . The question of parity for the 
remaining cases m = 4,8,12,... is a subject of current 
investigation. 



ACKNOWLEDGMENTS 

VVA thanks Hoda Hossein-Nejad for fruitful discus- 
sions. This work was supported by the Fulbright Canada 
Program and by the U.S. Air Force Office of Scientific 
Research under contract number FA9550-1 0-1-0260. 



[1] R. L. Fulton and M. Gouterman, 

|J. Chem. Phys. 35, 1059 (1961)| 
[2] E. T. Jaynes and F. W. Cummings, 

|P. IEEE 51, 89 (1963)1 
[3] Q.-H. Chen, T. Liu, Y.- Y. Zhang, and K.-L. Wang, 

ArXiv larXiv:1011.3280vll 
[4] B. W. Shore and 
| J. Mod. Optic. 40, 1195 (1993)| 



Phys. Scripta 76, 146 (2007) ; Y 
R Yu, Q. Liang, J.-QT Liang 



Phys. Rev. A 83, 065802 (2011) 



L. Knight, 
J. Larson, 
Zhang, G. Chen, 
and S. Jia, 
M. H. Naderi, 



J. Phys. A: Math. Theor. 44, 055304 (2011)1 



[5] T. Liu, K.-L. Wang, and M . Feng, 

|EPL (Europhysics Letters) 86, 54003 (2009) 
[6] B. Gardas, J. Phys. A: Math. Theor. 44, 195301 (2011). 

[7] V. V. Stepanov G. Miiller, and 

J. Stolze, |Phys. Rev. E 77, 066202 (2008) | 

P. L. Hagelstein and L U. Chaudhary, 
|J Phys. B: At. Mol. Opt. Phys. 41, 035601 (2008)] 
[8] C~ Emary and _R. F. Bishop, 

A. Pereverzev and 



J. Math. Phys. 43, 3916 (2002) 



E. R. Bittner,|Phys. Chem. Chem. Phys. 8, 1378 (2006) 



T. Werlang, A. Dodonov, E. Duzzioni, and C. Villas- 



Boas, |Phys. Rev. A 78, 053805 (2008)| I. Travenec and 
L. Samai. IArXivlarXiv:1107.4479l 
[9] U. Herfort and M. Wagner, J. Phys.: Condens. Matter 
13, 3297 (2001). 
[10] A. Eisfeld, L. Braun, W. T. Strunz, J. S. Briggs, J. Beck, 



and V. Engel, |J. Chem. Phys. 12 2, 134103 (2005) 
[11] C. J. Gan and H. Zhen g ,|Eur. Phys. J. D 59, 473 (2010) 



[12] S. Ashhab and F. Nori, |Phys. RevT A 81, 042311 (2010; 
[13] P. Forn-Dfaz, J. Lisenfeld, Eh Marcos, T! Garcia 
Ripoll, E. Solano, C. Harmans, and J. Mooij 
|Phys. Rev Lett. 105, 237001 (20lo7 



[14] A.-S. F. Obada, H. Hessian, and A.-B. A. Mohamed 



Physica E 42, 12 62 (2010) ; J. Hausinger and M. Grifoni, 
Phys. Rev. A 83, 030301 (R) (2011)] 

and H. R. Jauslin, 



[15] M. Amniat-Talab, S. Guerin, 
|J. Math. Phys. 46, 042311 (20057 



[16] I. D. Feranchuk, I. I. Komarov, and A. P. Ulyanekov, J. 
Phys. A: Math. Gen. 29, 4035 (1996). 



[17] E. Irish, | Phys. Rev. Lett. 99, 173601 (2007) 
[18] C. Ciuti 



Bastard, 



and I. Carusotto, 



|Phys. Rey. B 72, 115303 (2005) 



voret, 



M. Girvin, 



M. De- 
and R. J. Schoelkopf, 

Ann. Phys. (Leipzig) 16, 767 (2007) | R. J 

Schoelkopf and S. M. Girvin, |Natiire 451, 664 (2008) | 
J. Bourassa, J. Gambetta, A. Abdumalikov, 
O. Astafiev, Y. Nakamura, and A. Blais, 
|Phys. Rev. A 80, 032109 (2009) | C. M. Wilson, 
G. Johansson, T. Duty, F. Persson, M. Sandberg, and 
P. Delsing, |Phys. Rev. B 81, 024520 (2010)[ D. J. Al- 
ton, N. P. Stern, T. Aoki, H. Lee, E. Ostby, K. J. Vahala, 



[19] 



and H. J. Kimble, [Nature Physics 7, 159 (2011) 
T. Niemczyk, F. Deppe, H. Huebl, E. P. Menzel, 
F. Hocke, M. J. Schwarz, J. Garcfa-Ripoll, D. Zueco, 
T. Hummer, E. Solano, A. Marx, and R. Gross, 



|Nature Physics 6, 772 (2010)] 

[20] S. Phoenix, |J. Mod. Optic. 36, 1163 (1989)| 
[21] J. Casanova, G. Romero, I. Lizuain, J. Garcfa-Ripoll, 

[22] 

[23] 



and E. Solano, |Phys. Rev. Lett. 105, 263603 (2010) 
J. 



Hausinger 



and 



|Phys. Rev. A 82, 062320 (2010) 
Z.-G 



Jing, 



Lii, 



|Phys. Rev. A 79, 044305 (2009) 



M. 




Grifoni, 


and 


Z. 


Ficek, 




D. 


Zueco, 


and 


P. 


Hanggi, 



G~! Reuther, ST Kohler 

Phys. Rev. A 80, 033 846 (200 9)] M.-J. Hw ang and 
M.-S. Choi, |Phys. Rev. A 82, 025802 (2010) | B. Per- 
opadre, P. Forn-Dfaz, E. Solano, and J. Garcfa-Ripoll, 



Phys. Rev. Lett. 105, 023601 (2010) | Y.-Y. Zhang, Q.- 
H. Chen, and S.-Y. Zhu. IArXivlarXiv:1106.2T9"ri 
F. Beaudoin, J. M. Gambetta, and A. Blais, 



IArXivlarXiv :1107.3990 
[24] H. Zheng, S. Zhu, 



and M. S. Zubairy, 



|Phys. Rev. Lett. 101, 200404 (2008)| 
[25] I. Lizuain, J. Casanova, J. Garcfa-Ripoll, J. Muga, and 

E. Solano, |Phy~s. Rev. A 81, 062131 (2010)] 
[26] Q. Ai, Y. Li, H. Zheng, and C. P. Sun, 



8 



[27 

[28 

[29 

[30 

[31 
[32 

[33 

[34 

[35 
[36 
[37 



|Phys Rev. A 81, 042116 (2010)| Q. Ai and J.-Q. 
Liao, |Commun. Theor. Phys. 54, 985 (2010)| X. Cao, 
J. You^ PL Zheng, ~K~. Kofman, and F. Nori, 
|Phys. Rev. A 82, 022119 (2010) 



[38 
[39 

[40 

[41 
[42 
[43 
[44 

[45 
[46 

I 47 
[48 
[49 
[50 



R. 



R. 



Puri 



and 



R. 



K. 



Bullough, 
T. Dung 



|J. Opt. Soc. Am. B 5, 2021 (19 88) H 
and N. D. Huyen, Phys. Rev. A 49, 473 (1994). 
A. H. Toor and M. S. Zubairy, Phys. Rev. A 45, 4951 
(1992). 

C. Emary and R. F. Bishop, 

|J Phys. A: Math. Gen. 35, 8231 (2002)] 
C. Emary, On the Quantum Mechanics of Generalised 
Rabi Systems, Ph.D. thesis, UMIST, Manchester (2001). 



S. N. Dolya, |J. Math. Ph 


ys. 49, 053524 (2008) 


|J. Math. Phys. 50, 03351272009) 





B. 

(1981). 
J. L. 



Skinner 



and 



J. Phys. Chem. 90, 4931 (1986) 



D. 

L. 



132 

Hsu, 
Skinner, 



Annu. Rev. Phys. Chem. 39, 463 (1988) 



A. M. Perelomov, Generalized Coherent States and Their 
Applications (Springer, Berlin, 1986). 
R. Sohn and M. S. Swanson, 

|J. Phys. A: Math. Gen. 38, 2511 (2005)| 
M. Brune, J. Raimond, P. Goy, L. Davidovich, 
and S. Haroch e, |Phys. Rev. Lett. 59, 1899 (1987) 
P. Machnikowski, |Phys. Status Solidi C 5~2 486 (2008) 
M. Ornigotti, G. D. Valle, T. T. Fernandez, 
A. Coppa, V. Foglietti, P. Laporta, and S. Longhi, 
|J. Phys. B: At. Mol. Opt. Phys. 41, 085402 (2008)] 
E. del Vailed S. Zippilli, F. R Laussy, 
A. Gonzalez- Tudela, G. Morigi, and C. Tejedor, 
|Phys. Rev. B 81, 035302 (2010)| Y. Ota, S. Iwamoto , 
N. Kumagai. and Y. Arakawa. IArXivlarXiv:1107.0372l 
E. Irish, J. Gea-Banacloche, I. Martin, and K. Schwab, 
|Phys. Rev. B 72 195410 (2005)| 

N. Narozhny, J. Sanchez-Mondragon, and J. H. Eberly, 
|Phys. Rev. A 23, 236 (198T)] M. O. Scully and M. S. 
Zubairy, Quantum Optics (Cambridge University Press, 
Cambridge, 1997). 

M. Wagner, Unitary Transformations in Solid State 
Physics (North- Holland, Amsterdam, 1986). 
A. V. Leonov and I. D. Feranchuk, J. Appl. Spectrosc. 
77, 899 (2011). 

P. Feng, Y. Y ou-Kun, X. Ming-Xia, H. Wen- Juan, and 
J. P. Draayer, [Commun. Theor. Phys. 48, 53 (2007)] 
C. F. Lo, K. L. Lui, K. M. Ng, and P. H. Yuen, Quantum 
Semicl. Opt. 10, L63 (1998). 

H. B. Shore and L. M. Sander, Phys. Rev. B 7, 4537 
(1973); M. Wagner, J. Phys. A: Math. Gen. 17, 2319 
(1984). 

S. Paganelli and S. Ciuchi, 



J. Phys.: Condens. Matter 18, 7669 (2006) 



F. Pan, X. Guan, Y. Wang, and J. P. D raayer , 
|J. Phys. B: At. Mol. Opt. Phys. 43, 175501 (2010)] 
R. Graham and M. Hohnerbach, Z. Phys. B: Cond. Mat. 
57, 233 (1984). 

R. F. Bishop and A. Vourdas, Phys. Rev. A 50, 4488 
(1994). 

R. W. Munn and R. Silbey, J. Phys. A: Math. Gen. 11, 
939 (1978). 

V. S. Popov and A. M. Perelomov, Sov. Phys. 
JETP-USSR 29, 738 (1969); M. V. Satya- 



narayana, |Phys. Rev. D 32, 400 (1985)| P. Marian 
IPhys. Rev. A 44, 3325 (1991)' 



[51] C. F. Lo, K. L. Liu, and K. M. Ng, EPL (Europhysics 
Letters) 42, 1 (1998). 



9 



0.5 1 1.5 2 2.5 

A/a. 

FIG. 1: Comparison of the first few RWA (long blue 
dashes) and S-RWA (short red dashes) energies at 
4cj = J. The two features of the exact spectrum present 
in the S-RWA and not in the RWA are convergence of 
energies of opposing parity at large coupling and the 
ability of the S-RWA to properly retain the ground state 
for large values of A. The black circle depicts where the 
RWA ground state is surpassed by an excited state. 



0.2 0.4 0.6 0.8 1.0 1.2 1.4 

(a) A/a, 




0.5 1.0 1.5 2.0 

(b) 

FIG. 2: (a) Numerical energies (solid black) compared 
with the 2-by-2 S-GRWA (short red dashes) and GRWA 
(long blue dashes) at the resonance case 2 J = u = 1, 
where the differences are most evident, (b) The 4-by-4 
S-GRWA matrix truncation (dashed) plotted against 
the numerical energies (solid). The initial physical 
reasoning for the respective RWA and S-RWA 
approximations suggest that the GRWA should be used 
for the small coupling regime while the S-GRWA, with 
its respect for parity and a corrected ground state, 
should be used for A > 0.5. 
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(a) Voi 




0.1 0.2 0.3 0.4 0.5 



(b) V" 

FIG. 3: (color online) Numerical energies (solid black) 
compared with the S-GRWA energies (dashed) (a) at 
the resonance case J = uj = 1 and (b) at 2 J = ui = 1. 
The even (odd) manifold energies are red (blue); the 
respective positive (negative) parities on each manifold 

arc in large (small) dashes. Although the numerical 
energies have the same parities and come from the same 

respective even/odd manifolds as the analytical 
approximation, all numerical results are in black for an 
easier visual comparison to the approximation. The 

numerical energies have been obtained by 
approximately diagonalizing the four matrices in Eq. 
l[2"gj) truncated to 100. 
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TABLES 





A < w 


A > w 


ui x 2J 


RWA - 


> GRWA 


w > 2J 


S-GRWA <- 


- S-RWA 



TABLE I: Sketch of the RWA and S-RWA and their 
respective extensions in the parameter space. The 
arrows represent the direction of extension upon 
generalization of each respective approximation. 
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